Accepted for publication in ApJ 

Preprint typeset using I^'T^]X style emulatcapj v. 04/20/08 



O 

(N 

o 
O 



< 
o 



> 
oo 

d 



SUBMILLIMETER OBSERVATIONS OF DENSE CLUMPS IN THE INFRARED DARK CLOUD G049.40-00.01 
Miju Kang^, Minho Choi^, John H. Bieging^, Jeonghee Rho^, Jeong-Eun Lee", and Chao-Wei Tsai^ 

^ Korea Astronomy and Space Science Institute, 776 Daedeokdaero, Yuseong, Daejeon 305-348, Republic of Korea; mjkang@kasi.re.kr. 
^ Steward Observatory, University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721, USA 
3 SOFIA Science Center, USRA/NASA Ames Research Center, Moffet Field, CA 94035, USA 
Department of Astronomy and Space Science, Kyung Hee University, Yongin, Gyeonggi 446-701, Republic of Korea 
^ Infrared Processing and Analysis Center, California Institute of Technology, Pasadena, CA 91125, USA 

Accepted for publication in ApJ 

ABSTRACT 

We obtained 350 and 850 fim continuum maps of the infrared dark cloud G049. 40-00. 01. Twenty- 
one dense clumps were identified within G049. 40-00. 01 based on the 350 /im continuum map with 
an angular resolution of about 9'.'6. We present submillimeter continuum maps and report physical 
properties of the clumps. The masses of clumps range from 50 to 600 Mq. About 70% of the clumps 
are associated with bright 24 /im emission sources, and they may contain protostars. The most massive 
two clumps show extended, enhanced 4.5 /im emission indicating vigorous star-forming activity. The 
clump size-mass distribution suggests that many of them are forming high mass stars. G049. 40-00. 01 
contains numerous objects in various evolutionary stages of star formation, from pre-protostellar 
clumps to H II regions. 

Subject headings: ISM: individual objects (G049. 40-00. 01) — ISM: structure — stars: formation 



1. INTRODUCTION 

Massive st ars form in co l d, den se molecular clouds and 
in clusters (jLada fc Ladal 120031) . Infrared dark clouds 
(IRDCs) are complexes of cold (T < 20 K), dense {n > 
10^ cm~^) molecular gas, some of which are believed 
to be the progenitors of massiv e stars and star clusters 
(lEgan et al.l 119981: iCarev et all 119981: iRathborne et al.l 
|2006[) . IRDCs were identified as dark extinction fea- 
tures because the cold dust in IRDCs absorb s the bright 
mid-infrared emission of the Galactic plane (jEgan et al.l 
[1998,) ■ Cold, dense molecular gas and dust in IRDCs 
were confirmed based on obs ervations at milhmeter 
and submillimeter w avelengths (jCarev et al.lll998l 120001 : 



IRathborne et al.|[20b 5. 200 6]) IRDCs fr agment into mul- 
ti£le cores ([Beuther &: Henningj [20091) . iBattersbv et al.l 
([2OIQ1) divided the evolutionary sequence of IRDCs into 
four stages from a quiescent clump to an embedded 
H II region by combining millimeter-centimeter contin- 
uum data and spectroscopic data of the HCO+ and 
N2II+ lines. Some IRDCs are thought to be good tar- 
gets for investigatin g the initial conditions of massive star 
formation (Beuthe r et al.l I2007t IRathborne et al.l 120061: 

iChambers c t al. 20091). 

Recently, iKang et al.l (|2009[ ) catalogued embedded 
young stellar objects (YSOs) near W51 using the data 
from the Galactic Legacy Infrared Mid-Plane Su rvey Ex- 
traordinaire (GLIMPSE LiBenj amin et al ."2003') and the 
Mid-infrared Imaging Ph otometer for Spit zer Galactic 
plane Survey fMIPSGAL: ICarev et al.l[2009l) . A total of 
35 YSOs have been found over the area of 0? 15 x 0? 14 cen- 
tered at I = 49?4 and b = 0?00. This region corresponds 
to MSXDC G049.40-00.01 identified from the MSX 8 fim 
data (Simon ct al. 2006a). (Hereafter we refer the region 
to G049.40-00.01, after dropping the MSXDC label.) 
G049. 40-00. 01 includes three Spitzer dark clouds cata- 
logue d with the GLIMPSE 8 fim data (jPeretto &: Fulleii 
I2009D . G049.40-00.01 is associated with the CO emission 
showing a velocity peak at Vlsr = 61 km s~^, which is a 



part of the "cl uster" region near the active star forming 
complex W51 (jKang et al.ll20ldl ). Recent measurements 
of the distance to W51 gave 5.4lto;28 ^P*^ (jSato et al.l 
[ Mg), 5.i;i;^ kp c (IXu et al.l [2001) . and 6.1 ± 1.3 kpc 
(jlmai et al.ll2002f ). Here, we adopt 6 kpc as the distance 
to G049. 40-00. 01 (with an uncertainty of ~1 kpc), for 
consistency with previous works. The peak H2 column 
density estimated from the ^ ^CO J = 1 — li ne ob- 
servations is 1.6 X 10 cm" -2 (iSimon et al.ll2006b[ ). Two 
comp act H II regions we re identified based on the Spitzer 
data ([Phillips &: Ramos-Larios 2008 ). 

In this paper, we present the results of submillimeter 
observations of the IRDC G049.40-00.01 using the Sub- 
millimeter High Angular Resolution Camera II (SHARC- 
II) at the Caltech Submillimeter Observatory (CSO). We 
describe details of the observations and data in Section 
2. Then we report the results in Section 3 and discuss 
the physical properties of the clumps in G49. 40-00. 01 in 
Section 4. We summarize the main results in Section 5. 

2. OBSERVATIONS AND DATA 

The observations were made on 2010 May 10 and 15 
at the CSO 10.4 m telescope near the summit of Mauna 
Ke a, Hawaii. We used the bolometer camera SHARC- 
II (jDowell et al.l 120031 ) . The instrument resolution of 
SHARC-II is 8'.'0 at 350 //m and 19'.'4 at 850 fim. The 
Dish Surface Optimization System was used to correct 
the dish surface for static imperfections and deforma- 
tions d ue to gravitational forces as the dish moves in ele- 
vation (jLeong et al.ll2006[ ). We obtained five scans at 350 
/im for a total integration time of 50 minutes in moder- 
ate weather (t225GHz ~ 0.046-0.068) and fourteen scans 
at 850 /im for a total integration time of 140 minutes 
in moderate weather (t225GHz ~ 0.053-0.078). Point- 
ing and calibration scans were taken on hourly basis on 
strong submillimeter sources: Neptune, Arp 220, IRAS 
16293-2422, CRL 2688, W75N, and K 3-50. The data 
were reduced using version 2.01-2 of the software package 
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Fig. 1. — Infrared and submillimcter maps of the IRDC G049. 40-00. 01. (a) IRAC three-color composite (8.0 fira in red, 4.5 /im in green, 
and 3.6 /im in blue) image overlaid with the Aa contour of the 350 /im map. Two large crosses mark the compact H II regions PR 29 
and 30, from north to south UPhillips fc Ramos-Lario3l200SI) . (b) IRAC and MIPS two-color composite (8.0 /im in cyan and 24 fim in 
red) image. The scale bar indicates 2 pc. (c) SHARC-II 350 fira image (both gray scale and contours). The compact sources (clumps) 
are labeled and marked by small crosses. Contour levels are 4, 8, 12, 16, and 20 X 0.09 (Ict) Jy beam~^. (d) SHARC-II 850 iim image. 
Contour levels are 3, 6, 9, and 12 X 0. 05 (Ict) Jy beam-^ The crosses mark the positions of the 350 fim clumps. The filled triangles 
mark the YSOs from [Kang et aLl II2009I V The FWHM resolutions of the SHARC-II maps are indicated in the bottom left corners. The 
extraneous contours (noisy areas) near the edges of the SHARC-II maps should be ignored. 



CRUSH (|Kovacsll2006l Comprehensive Reduction Utility 
for SHARC-II). The final maps were smoothed to angu- 
lar resolutions of FWHM = 9'.'61 at 350 ^m and 23'.'32 
at 850 /im. 

The Spitzer data presented here are the GLIMPSE 
and MIPSGAL data. The GLIMPSE and MIPSGAL 
are legacy programs covering the inner Galactic plane at 
3.6, 4.5, 5. 8, and 8.0 iim wit h the Infrared Array Cam- 
era (IRAC: lFazio et al.|[200l ) and at 24 and 70 ^m with 
the Multiband Im aging Photometer for Spitzer (MIPS; 
iRieke et"al]l2004[ ). respectively. 

3. RESULTS 

Figure [T] shows the maps of the Spitzer IRAC, MIPS, 
and SHARC-II 350 and 850 ^im toward the IRDC 
G049. 40-00. 01. In the IRAC composite map (Figure 
W^), a dark filamentary structure is seen in absorption 
against the background emission. The bright features in 
the IRAC composite map are two H II regions and bright 
diffuse emission from the active star- forming region W51 
to the south. The spatial distribution of the dark fil- 
amentary structure agrees well with the distribution of 



the 350 /Ltm emission. Figure lb shows the two-color com- 
posite (IRAC 8.0 in cyan and MIPS 24 fim in red) 
image, overlaid with the 350 fim emission map. Some 
regions of G049. 04-00. 01 are stiU dark at 24 fim, imply- 
ing low temperature and high column density. There are 
also many bright 24 ^m point-like sources located within 
the IRDC. The dark features of the IRDC seen in the 
Spitzer maps appear in emission in the 350 and 850 /im 
maps (Figures Ic and Id) because their cold thermal dust 
emission peaks at millimeter/submillimeter wavelengths. 
The bright clumps in the 350 /xm map are well aligned 
with the 24 fim peaks, which suggests that the 24 fim 
sources are central stars of these dense clumps. The dis- 
tribution of 850 /im emission is very similar to that of 
350 /im emission except the details unresolved with the 
relatively large beam of the 850 fim map. The maximum 
flux densities are 1.96 ± 0.09 and 0.64 ± 0.05 Jy beam^i 
at 350 and 850 fim, respectively. 

A total of 21 sources were identified as compact clumps 
for which peak flux densities are greater than the Aa level 
in the 350 fim continuum map by eye. When a peak is 
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TABLE 1 

SUBMILLIMETER CONTINUUM SOURCE PARAMETERS 



Peak Position^ 
I b 



Peak Flux'' 



850 Atm 



350 fiia 



Total Flux 
at 350 fim 



Size" 



24 fim 
Sources'* 



Associated 
YSOs<= 



Source 


(.deg) 


(deg) 


(Jy beam ^) 


(Jy beam ^) 


(Jy) 


( 


X 


] 






1 


49.361 


0.037 


0.30 


0.58 


3.4 ± 0.4 


21 


X 


21 






2 


49.366 


0.005 


0.24 


0.54 


1.6 ± 0.3 


21 


X 


8 


2 


345, 348 


3 


49.367 


0.025 


0.58 


1.75 


11.8 ± 1.1 


22 


X 


19 


2 


346 


4 


49.390 


-0.016 


0.26 


0.54 


2.5 ± 0.2 


26 


X 


21 


3 


362, 366, 367 


5 


49.393 


0.019 


0.23 


0.55 


5.1 ± 1.6 


53 


X 


19 






6 


49.396 


-0.000 


0.23 


0.41 


1.2 ± 0.1 


16 


X 


15 


1 




7 


49.401 


-0.036 




0.54 


1.6 ± 0.2 


22 


X 


16 






o 

O 


/in /I HQ 


U.UUo 


V.Zi 


U.OD 


4.4 ± U.D 


A Q 
4o 


X 


OQ 
ZO 






9 


49.408 


-0.038 


0.19 


0.68 


2.3 ± 0.3 


29 


X 


13 


1 




10 


49.409 


-0.007 


0.48 


1.09 


5.5 ± 1.1 


21 


X 


18 


3 




11 


49.410 


-0.016 




0.63 


1.4 ± 0.2 


16 


X 


11 






12 


49.412 


-0.020 




0.99 


2.2 ± 0.1 


20 


X 


11 


1 




13 


49.416 


-0.010 




0.57 


1.4 ± 0.3 


21 


X 


11 


1 




14 


49.418 


-0.038 




0.88 


3.4 ± 0.3 


29 


X 


13 


1 




15 


49.417 


-0.027 




0.70 


2.1 ± 0.5 


14 


X 


13 






16 


49.420 


-0.021 


0.65 


1.96 


13.3 ± 0.8 


27 


X 


20 


3 


409, 414, 419 


17 


49.423 


-0.040 


0.41 


1.44 


5.6 ± 0.8 


24 


X 


17 


1 


421 


18 


49.425 


-0.013 


0.27 


1.02 


4.3 ± 0.3 


33 


X 


16 


1 f 




19 


49.434 


-0.034 


0.21 


0.72 


1.2 ± 0.3 


13 


X 


7 






20 


49.439 


-0.038 


0.26 


1.31 


3.4 ± 0.9 


14 


X 


12 


1 t 




21 


49.443 


-0.032 


0.40 


0.94 


2.3 ± 0.4 


21 


X 


17 


1 


438 



^ Peak position in tiie 350 /im image. 

^ Uncertainties are 0.05 and 0.09 Jy bcani"^ at 850 and 350 /ini, respectively. 

" Deconvolved FWHM size of each clump in the 350 ^m map. The first number is the size along the major axis (longest 
diameter), and the second number is the size along the direction perpendicular to the major axis. 
^ Number of 24 ^im point sources within the size of each clump. 
" YSO number in Table 2 of IKang et~all ||2009(1 . 



Associated with an H II region. 

located close to another one, it is considered an inde- 
pendent peak if the peak intensity is higher than the 
intensity at the interface between them by Icr level or 
larger and if their separation is larger than a beam size. 
Table [T] lists the peak position, peak flux densities, 350 
/zm total flux density, size, number of associated 24 ^m 
sources, and associated YSOs of each clump. The total 
flux density at 350 /im is measured from the emission 
in the circumscribed box of 2a contour. In crowded ar- 
eas, the saddle point between adjacent clumps is used 
to limit the area for total flux measurements. The size 
represents the FWHM, deconvolved with the beam. In 
crowded areas, the intensity at the saddle point between 
adjacent clumps can be higher than the half- maximum 
intensity. In this case, we give the distance to the saddle 
point from the peak position as a size (in the direction 
to the adjacent clump). Column 8 lists the number of 
all 24 i^m point sources detected with significant signal- 
to-noise rati os (S/ N > 4). Column 9 lists the YSOs in 
|Kang et al] (j2009t ). They used a detection criterion of 
S/N > 7 at 24 fim to identify YSOs. The peak fluxes at 
850 /im are measured only for the clumps having 850 fxin 
peaks within the size boundaries of the 350 /im sources 
as described above. Table [2] lists the 850 /im total flux 
densities. For many clumps the 850 /im total flux cannot 
be measured individually because of crowding. 

4. DISCUSSION 

4.1. Physical Parameters 

Continuum emission at 350 /im is sensitive to dust 
temperature, and the mass of molecular gas estimated 
from the 350 flux density depends on the assumed 
dust temperature that varies with the presence or ab- 



TABLE 2 
Total flux at 850 fim 



Sources 


Total Flux 

(Jy) 


1 


0.30 ± 0.05 
0.29 ± 0.13 
1.03 ± 0.15 
0.26 ± 0.05 


2 


3 


4 


5, 6, 8 


0.98 ± 0.18 


7, 9 


0.22 ± 0.12 


10, 13 


0.79 ± 0.15 


11, 12, 15, 16, 18 


1.44 ± 0.18 


14, 17 


0.73 ± 0.15 


19, 20, 21 


0.51 ± 0.14 


Note. — For crowded 
flux densities of each area 


areas, the total 
are listed. 



sence of central heating by embedded protostars. The 
clumps without detectable 24 /im so urces may be pre- 
protos tellar (harboring no protostar). iHennemann et ahl 
(|2009D derived dust temperatures of 22 and 1 5 K for cores 
with a nd without 24 /im source, respectively. iStutz et al.l 
([2010) found that the dust temperatures are ^^17. 7 K 
near the protostar and ^10.6 K for the starless core in 
the Bok globule CB 244. Pre-protostellar and proto- 
stellar cores in the IRDC GOll. 1-0.12 have core tem- 
peratu res of 22 K (jHenning et al.ll2010[ ). iWilcock et all 
(|2011f) derived temperatures of 8-11 K at the center 
of the cores and 18-28 K at the surface using radia- 
tive t ransfer models of IR DC seen in Herschel observa- 
tions. iPeretto et al.l ()201ClD reported that IRDCs are not 
isothermal, showing that the dust temperature decreases 
significantly within IRDCs, from background tempera- 
tures of 20-30 K to minimum temperatures of 8-15 K 
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TABLE 3 

Properties of the Clumps 



Source 



Mtotal" 



MpwHM 



^ Concentration 



1 


154 


± 


18 


149 


± 


27 


0.65 


It 


0.18 


2 


73 


± 


14 


54 


± 


20 


0.73 


it 


0.44 


3 


536 


1 

± 


50 


354 


± 


75 


0.77 


1 

± 


0.05 


4 


113 


1 


9 


127 


± 


14 


0.65 


1 

± 


0.21 


5 


232 


± 


73 


226 


± 109 


0.78 


it 


0.54 


6 


54 


1 

± 


5 


58 


± 


7 


0.57 


1 

± 


0.17 


7 


73 


± 


9 


88 


± 


14 


0.50 


it 


0.08 


8 


200 


± 


27 


217 


± 


41 


0.55 


it 


0.09 


9 


104 


± 


14 


99 


± 


20 


0.59 


it 


0.08 


10. . . . 


250 


± 


50 


182 


± 


75 


0.72 


it 


0.13 


11 ... . 


64 


± 


9 


60 


± 


14 


0.50 


it 


0.07 


12.... 


100 


± 


5 


103 


± 


7 


0.64 


it 


0.17 


13. . . . 


64 


± 


14 


64 


± 


20 


0.39 


it 


0.06 


14. . . . 


154 


± 


14 


171 


± 


20 


0.69 


it 


0.23 


15.... 


95 


± 


23 


90 


± 


34 


0.62 


it 


0.13 


16. . . . 


604 


± 


36 


455 


± 


54 


0.67 


it 


0.03 


17. . . . 


254 


± 


36 


242 


± 


54 


0.71 


it 


0.14 


18. . . . 


195 


± 


14 


231 


± 


20 


0.54 


it 


0.05 


19. . . . 


54 


± 


14 


61 


± 


20 


0.62 


it 


0.08 


20 ... . 


154 


± 


41 


139 


± 


61 


0.64 


it 


0.04 


21 ... . 


104 


± 


18 


131 


± 


27 


0.65 


it 


0.08 


Note. 


— The mass 


is based on the total f 


ux 


density 



at 350 /xm assuming — 15 K. If T^i — 10 K, multiply 
by 4.1. If Td = 20 K, multiply by 0.47. 
^ Mass within the 2cr contour. 

^ Mass within the FWHM boundary with the opaeity 
scaled by 1/1.5. 

within the clouds. In this paper, we assume a dust tem- 
pera ture of 15 K for simplicity and for ease of comparison 
with lKauffmann fc P iUai (2010) (see Section 4.3). 

The masses of the submillimeter clumps in G049.40- 
00.01 were calculated using the 350 //m flux densities. 
The mass can be estimated by M — F^D'^ / K^Bj^iTd), 
where F^, is the flux density at 350 fim, D is the dis- 
tance to the source, B^, is the Planck function, Td is the 
dust temperature, and is the dust mass opacity. The 
value of at 350 yum used to calculate the masses is 
5.91 cm^ from the coagulated dust m odel with thin 
ice mantles of lOssenkopf fc Henninj (|1994[ ). Table [3] lists 
the total mass within the 2a contour. The representative 
masses listed in Table [3] were calculated using a uniform 
temperature of 15 K. These masses can be underesti- 
mates for colder clumps (those without a 24 /im source) 
and overestimates for hotter clumps (those with 24 /im 
sources). Different assumptions on the dust temperature 
can increase or decrease the mass estimate by a factor of 
~3 (see the footnote of Table [3l). The uncertainty of the 
mass caused by the uncertainty of the distance is ~30%. 
The sum of the clump masses is ^3,600 M©. 

4.2. Comparison with Simple Models 

The clumps in G049. 40-00. 01 could be produced by 
gravitational fragmentation, and some simple quantities 
can be calculated to check the consistency. The average 
surface density is ~0.04 g cm~^ for the area within the Icr 
contour in the 350 /im map. With this surface density 
and the assumed dust temperature of 15 K, the criti- 
cal wavelength of an isothermal, infinite, self-gravitating 
cylinder is Ac « 0.3 pc , and the corresponding critical 
mass is Mc « 5 Mq (jHartmannl l2002t iLarsonl [l985l ) . 
For comparison, the projected distance between the near- 
est neighbours of clumps in G049. 40-00. 01 lies between 
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Fig. 2. — Mass-concentration diagram for the clumps detected 
in the 350 fim map. Clumps with large (> 0.2) uncertainties 
in concentration are omitted. Red and green circles represent 
clumps associated with 24 fim sources and green fuzzies, respec- 
tively (see Section 4.3). The vertical dotted lines (C = 0.33 and 
0.72) represent the lim iting cases of stable Bonnor-Ebert spheres 
II Johnstone et al.|[2000D . 



0.4 and 2 pc. The mean clump separation is 0.9 pc, 
and the average clump mass is 170 M©. Many filamen- 
tary IRDCs fra gment into clumps with a similar frag- 
ment ation scale (jHenning et al.|[2010HMiettinen fc Hariul 
|2010[) . The difference between the calculated critical 
mass and the measured average mass probably indicates 
that the fragmentation history of G049. 40-00. 01 is more 
complicated than the fragmentation of an idealized sim- 
ple cylinder. Previously suggested possibilities include an 
extra support by turbulence, changes of cl oud tempera- 
ture, and a compression by ext ernal forces (|Onishi et al.l 
[1991 iMiettinen fc Hariul[20lol) . 

To study further physical conditions of the clumps, 
we compare the clumps to a static cloud model. 
The Bonnor-Ebert model describes the simplest self- 
gravitating pressure-c onfined isoth e rmal sphere in a hy- 
drostatic equilibrium (jEbertl 119551: iBonnod [19561) . The 
degree of self-gravitation within each clump can be esti- 
mated by calculating the degree of concentration of each 
clump. The concentration can be defined by C = l — Rs, 
where Rs is the ratio of average to central colu mn den- 
sity (see equation (4) of iJohnstone et all ()2000D ). Small 
concentrations (C < 0.33) imply uniform-density non- 
self-gravitating objects, and large val ues (C > 0.72) 
imply critically self-gravitating objects (jJohnstone et al.l 
l2000t ). The concentrations of the submillimeter clumps 
are listed in Table [H Figure [2] shows the concentration 
of the clumps against the clump mass. Most clumps 
in G049. 40-00. 01 have concentrations between 0.33 and 
0.72. Clump 3 (C = 0.77 ± 0.05) is more concentrated 
than what is permitted by stable Bonnor-Ebert spheres. 
In general, massive clumps seem to have high concen- 
trations. It is interesting to note that the four highest- 
concentration clumps are also the highest-mass ones, and 
they are all associated with 24 /im sources, indicating 
ongoing star formation. Two of them also show signs of 
shocked gas (see Section 4.3). For the rest of the clumps. 
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Fig. 3. — Mass-size diagrams for the clumps detected in the 
350 fira map. (a) Mass-size diagram with the total mass given 
in Table [3] (b) Mass-size diagram with the mass estimated from 
the flux density within the FWHM boundary (and with t he dust 
ogacity scaling factor applied, as suggested by Kauffmann fc Pillail 
12010 ')'). Clumps associated with 24 fim sources arc marked by red 
circles. Two large green circles mark the "green fuzzies" . The solid 
line represents the mass-size threshold fo r ma ssive star formation, 
mKp, proposed bv lKauffmann fc Pillail f2010l) . 

there is no clear difference between those with and with- 
out 24 /xm sources. 

4.3. High-mass Star Formation 

Massive stars are expected to be formed in molec- 
ular clouds with a large mass concentrated in a rela- 
tively small volume. Kauffmann & Pillai (2010) sug- 
gested a threshold for massive star formation by com- 
paring clouds with and without massive star formation: 
TOkp = 870 Mq r^-^^, where r is the effective radius (half 
of the geometric mean of the FWHM sizes in Table 1) 
in pc. Clouds/clumps more massive than this thresh - 
old seem to form massive stars. iParmentier et al.l ()2011l ) 
provided an explanation for the KaufFmann-Pillai thresh- 
old by calculating the probable mass of the most massive 
star formed in model clumps having power-law density 
profiles. 

Figure [3] shows the mass -size relation for the clumps 
in G049.40-00.01. KauffmaHEI&^illail (2010) assumed 
a dust temperature of 15 K. They used a dust opacity 
scaled by a factor of 1/1.5 (to match the masses esti- 
mated from dust emission and extinction) and a mass 
correction factor of l/ln2 (to convert the total mass to 
the mass contained in the FWHM), and these factors 
cancel out almost exactly. Figure [3^ shows the relation 
between the total mass and the effective radius. Figure 
[3}3 shows the relation between the mass contained within 
the FWHM boundary (with the op acity scaled by 1/1.5 
to ma ke a proper comparison with iKaufFmann fc Pillail 
(|2010( )) and the effective radius. The masses are listed 
in Table [3l All the clumps are distributed near tokp, 
within a factor of '-~^3.5. 

Association with bright 24 fim point-like sources is an 
indicator of star formation because 24 /im emission traces 
warm dust heated by the material accreting onto a cen- 
tral protostar. Fourteen clumps out of 21 have 24 /Ltm 



point-like sources within their extent. The most massive 
ones (clumps 3 and 16) even show extended, enhanced 4.5 
/xm emission called "green fuzzies" or extended green ob- 
jects, which indicate sh ocked gas ([Chambers et al.ir2009l : 
iCyganows ki et aD l2008'). Extended 4.5 fim emission also 
indicates vigorous star-forming activity. In the G049.40- 
00.01 region, two compact H II regions (PR 29 and 
30 in Figure [1^) were identifie d based on Spitzer data 
()Phillips fc Ramos-Lario32008l) . PR 29 is more extended 
than PR 30 in the infrared maps (Figures la and lb) and 
is not directly associated with a submillimeter clump. 
PR 29 seems to be a blister type H II region formed on the 
cloud surface near clump 18. In contrast, PR 30 is closely 
associated with the submillimeter clump 20, which sug- 
gests that it is still embedded in the dense cloud. PR 30 
may be less evolved than PR 29. Clumps 3, 16, and 20 
(clumps associated with "green fuzzies" or compact H II 
region) are located well above mxp, and this fact seems 
to corrobora te the threshold for massiv e star formation 
suggested bv lKauffmann fc Pillail (|2010D . 

Several clumps are dark at 24 /im, and they may be 
clumps either in the pre-protostellar phase or containing 
protostars of undetectably low luminosity. These clumps 
reside below tukp in Figure [3] and may not be dense 
enough to form stars (yet). Therefore, objects in vari- 
ous stages of star formation (from pre-protostellar cores, 
protostars, to H II regions) are distributed in and around 
the IRDC G049.40-00.01. 

Spectral observations would be required to definitively 
determine the evolutionary status of the clumps. Re- 
cent multi- wavelength, high angular resolution studies on 
IRDCs focused on the chemical, kinematic, and physical 
properti es of the ini t ial co nditions for massive star for- 
mation. iPillai et al.l (|2011l ) investigated secondary mas- 
sive cold cores in the vicinity of H II regions based on the 
line (NH2D, NH3, and HCO"'') and millimeter continuum 
observations. They showed that the cores in the earli- 
est stage of massive star formation are col d, dense, and 
highly deuterated ([NH2D/NH3] > 6%). iDevine et all 
([2011) found an anti-correlation between NH3 and CCS 
toward the IRDC G19.30+0.07 based on interferomet- 
ric observations at 22 GHz. The different evolutionary 
states of the young objects in G049. 40-00. 01 can be in- 
vestigated in detail by obtaining high-resolution interfer- 
ometric data in the future. 

5. SUMMARY 

We observed the IRDC G049.04-00.01 in the 350 and 
850 /im continuum with the SHARC-II bolometer cam- 
era. The dark features in the infrared images have a 
good agreement with the emission structure in the sub- 
millimeter images. Twenty-one clumps were identified 
based on the 350 /im continuum map, and the mass of 
each clump ranges from 50 to 600 Mq. The majority of 
these clumps are associated with bright 24 fim emission 
sources, indicating star- forming activity. The most mas- 
sive clumps (clumps 3 and 16) show extended, enhanced 
emission in the IRAC 4.5 fim image. All the clumps 
are distributed near the threshold for massive star for- 
mation. The IRDC G049.04-00.01 contains objects in 
various evolutionary stages of star formation. 
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